CLAIMS 

What is claimed is: 



1 1 . A method of modeling failure of a borehole in a^subsurface formation, the 

2 method comprising:; 

3 (a) defining a subsurface model including a plurality of regions, said 

4 plurality of regions including the borehole and at least one additional 

5 region selected from (i) a liner in the borehole, (ii) a casing in the 

6 borehole, and (iii) at least one earth formation , each of said plurality 

7 of regions comprising a plurality of nodes interconnected by a 

8 plurality of linkages, 

9 (b) defining material properties associated with said nodes and said 

I o linkages of said subsurface model, said material properties having a 

I I statistical variation; 

12 (c) specifying an initial deformation pattern of the model; and 

1 3 (d) using a dynamic range relaxation algorithm (DRRA) to find a force 

14 equilibrium solution for said subsurface model and said initial 

1 5 deformation pattern giving a resulting deformed model including 

16 fracturing. 

1 2. The method of claim 1 , wherein said nodes are arranged in a grid that is one of 

2 (i) a triangular grid, and, (ii) a random grid. 

1 3. The method of claim 1 wherein said linkages are selected from the group 

2 consisting of (A) springs, (B) beams, and. (C) rods. 



CON- 1021 



18 



The method of claim 1 wherein said linkages comprise springs, the method 
further comprising defining a normal force associated with each spring. 



1 5. The method of claim 1 wherein said linkages comprise beams, the method 

2 further comprising defining at least one of (A) a normal force, and (B) a shear 

3 force associated with each beam. 

1 6. The method of claim 1 wherein said linkages comprise rods, the method 

2 further comprising defining at least one of (A) a normal force and (B) a force 

3 associated with an angle between pairs of said adjacent ones of the plurality of 

4 rods. 

1 7. The method of claim 1 , wherein using the dynamic range relaxation algorithm 

2 further comprises applying said initial deformation model in a plurality of 

3 steps, each step comprising applying a specified fraction of the initial 

4 deformation and determining if any linkages between the nodes have been 

5 deformed beyond a breaking point and identifying a subset of the linkages that 

6 have been so deformed. 

1 8. The method of claim 7, wherein applying the dynamic range relaxation 

2 algorithm further comprises iteratively breaking the one linkage of the subset 

3 of linkages that has been deformed the most and applying a relaxation 

4 algorithm to the remaining unbroken linkages. 
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The method of claim 9 wherein the at least one earth formation further 
comprise a near earth formation including a gravel pack and a far earth 
formation. 



1 10. The method of claim 1 wherein the plurality of regions comprises a liner in the 

2 borehole, an earth formation including a near earth formation and a far earth 

3 formation, and a gravel pack disposed between the liner and the near earth 

4 formation. 

1 11. The method of claim 1 wherein said linkages connect at least one selected 

2 node of said plurality of nodes with (i) a plurality of nearest neighbors of the 

3 at least one selected node, and (ii) a plurality of next nearest neighbors of the 

4 at least one selected node. 

1 12. The method of claim 1 wherein said earth formations include a fluid, said fluid 
: flowing into the borehole, and said deformation pattern is determined in part 

3 by a decrease in formation fluid pressure resulting from flow of said fluid into 

4 the borehole. 

1 13. The method of claim 1 2 wherein using the DRRA further comprises 

2 determining an additional force at each node related to a difference in said 

3 fluid pressure on opposite sides of at least a subset of the plurality of nodes. 
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1 14. The method of claim 13 wherein determining said additional force further 

2 comprises performing a simulation selected from (i) a finite difference 

3 simulation, and, (ii) a finite element simulation, of said fluid flow. 

1 15. The method of claim 14 wherein performing said simulation further comprises 

2 changing at least one of (A) a permeability, and, (B) a porosity used in said 

3 simulation responsive to said deformation. 

1 1 6. The method of claim 1 wherein said borehole includes a substantially vertical 

; " 2 section wherein said initial deformation pattern is substantially azimuthally 

; : j 3 symmetric about an axis of the borehole in said section. 

n 

1 17. The method of claim 16 wherein said borehole includes a deviated section 

* 2 wherein said initial deformation pattern is asymmetrical about an axis of the 

; u 3 borehole. 

1 18. A method of modeling failure of a borehole in a subsurface formation, the 

2 method comprising: 

3 (a) defining a subsurface model having a plurality of nodes and including 

4 a plurality of regions, said plurality of regions including the borehole 

5 and at least one additional region selected from (i) a liner in the 

6 borehole, (ii) a casing in the borehole, and (iii) at least one earth 

7 formation, each of said plurality of regions comprising a plurality of 

8 nodes interconnected by a plurality of linkages, 
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9 (b) defining material properties associated with said nodes and said 

10 linkages of said subsurface model, said material properties having a 

1 1 statistical variation; 

12 (c) specifying a force distribution applied to the model at boundary nodes 

1 3 of said plurality of nodes; and 

1 4 (e) using a dynamic range relaxation algorithm (DRRA) to find a force 

1 5 equilibrium solution for said subsurface model and said force 

1 6 distribution giving a resulting deformed model including fracturing. 

1 1 9. The method of claim 1 8 wherein the subsurface formation has been subjected 

2 to large scale geologic deformation and wherein specifying said force 

3 distribution further comprises: 

4 (i) simulating the large scale geologic deformation to determine a stress 

5 distribution in the subsurface formation in the absence of the borehole, 

6 (ii) defining a trajectory for the borehole therein, and 

7 (iii) identifying locations along said trajectory that are likely to fail. 

1 20. The method of claim 1 8 wherein the forces can vary between the boundary 

2 nodes. 

1 21. The method of claim 1 9 wherein identifying said trajectories further comprises 

2 removing a plurality of nodes along said trajectory. 

1 22. The method of claim 1 8, wherein said nodes are arranged in a grid that is one 
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of (i) a triangular grid, and, (ii) a random grid. 

The method of claim 1 8 wherein said linkages are selected from the group 
consisting of (A) springs, (B) beams, and. (C) rods. 

The method of claim 1 8 wherein said linkages comprise springs, the method 
further comprising defining a normal force associated with each spring. 

The method of claim 18 wherein said linkages comprise beams, the method 
further comprising defining at least one of (A) a normal force, and (B) a shear 
force associated with each beam. 

The method of claim 18 wherein said linkages comprise rods, the method 
further comprising defining at least one of (A) a normal force and (B) a force 
associated with an angle between pairs of said adjacent ones of the plurality of 
rods. 

The method of claim 1 8, wherein using the dynamic range relaxation 
algorithm further comprises applying said force distribution in a plurality of 
steps, each step comprising applying a specified fraction of the force and 
determining if any linkages between the nodes have been deformed beyond a 
breaking point and identifying a subset of the linkages that have been so 
deformed. 
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1 28. The method of claim 27, wherein applying the dynamic range relaxation 

2 algorithm further comprises iteratively breaking the one linkage of the subset 

3 of linkages that has been deformed the most and applying a relaxation 

4 algorithm to the remaining unbroken linkages. 

1 29. A method of modeling faulting and fracturing in a subsurface volume of the 

2 earth comprising: 

3 (a) defining said subsurface model including a plurality of interconnected 

4 nodes and material rock properties within the subsurface volume; 

5 (b) specifying a stress distribution at a subset of said plurality of nodes, 

6 said subset comprising boundary nodes; and 

7 (c) using a dynamic range relaxation algorithm to find a force equilibrium 

8 solution for said subsurface model and said stress distribution giving a 

9 resulting deformed model including fracturing. 

1 30. The method of claim 29, wherein defining a subsurface model, and specifying 

2 said stress distribution further comprises using a graphical user interface. 

1 31. The method of claim 29, wherein said nodes are arranged in a grid that is one 

2 of (i) a triangular grid, and, (ii) a random grid. 

1 32. The method of claim 29, wherein said nodes are interconnected by linkages 
selected from (i) springs, (ii) beams, and, (iii) rods. 
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